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Abstract. We review various combinatorial problems with underlying classical or quan- 
tum integrable structures. (Plenary talk given at the International Congress of Mathe- 
matical Physics, Aalborg, Denmark, August 10, 2012.) 



1. Introduction 

^ . In these notes, we present a few mathematical problems or constructs, most of them 

^ combinatorial in nature, that either were introduced to explicitly solve or better understand 

physical questions (Mathematical Physics) or can be better understood in the light of 
physical interpretations (Physical Mathematics). The frontier between the two is subtle, 
and we will try to make this more concrete in the examples chosen below. 

We start from two different physical theories of discretized random surfaces, describing 
the possible fluctuations of an underlying two-dimensional discrete space time or equiv- 
alently of a discretized Lorentzianjl] (Sect. 2) or Euclidian[38] (Sect. 3) 1-1-1- or 2- 
dimensional metric. In both simple models, we show the existence of some hidden integrable 
structure, in two very different forms. The first model may be equipped with an infinite 
family of commuting transfer matrices p5] governing the time-evolution of the surfaces, thus 
displaying quantum integrability, with an infinite number of conserved quantities. The 
second model actually addresses the correlations of marked points on the discrete surface 
at prescribed geodesic distance, and we show that such correlations viewed as evolutions 
in the geodesic distance variable, form a discrete classical integrable system with some 
conserved quantities modulo the equations of motion[8j. Remarkably, both cases can be 
reduced bijectively to statistical ensembles of trees. 

Few combinatorial objects form such a vivid crossroads between various physical and 
mathematical theories and constructs as the Alternating Sign Matrices[6]. These were 
defined by Robbins and Rumsey[35j in an attempt to generalize the notion of determinant 
of a square matrix, while keeping certain properties, in particular the Laurent polynomiality 
of the result in terms of the matrix elements. Interestingly, this happened two decades 
before a combinatorial theory of the Laurent phenomenon was discovered by Fomin and 
Zelevinsky, under the name of Cluster Algebra[2n]- In Sect. 4, we try to unravel the thread 
from alternating sign matrices to integrable lattice models, back to more combinatorial 
objects such as plane partitions. We show that the refined enumeration of such objects 
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involves underlying integrable structures similar to that of discrete Lorentzian gravity. This 
fact is actually used to prove the Mills Robbins Rumsey conjecture[53] relating alternating 
sign matrices to descending plane partitions. Section 5 is devoted to the cluster algebra 
formulation and explicit solution of the celebrated T-system|28j. a discrete 2+1-dimensional 
integrable system |27j related to alternating sign matrices, but also to domino tilings of the 
Aztec diamond [T8| |37] and Littlewood- Richardson coefficients[2l]. We show that solving 
such a system amounts to computing partition functions of weighted path models, a kind 
of discrete path integral, that gives a new insight on the Laurent positivity conjecture of 
cluster algebra. 



2. 1 + lD Lorentzian gravity, integrability and trees 

2.1. 1+lD Lorentzian gravity. Discrete models for I+ID Lorentzian gravity are defined 
as follows. They involve a statistical ensemble of discrete space-times, which may be 
modeled by random triangulations with a regular time direction (a segment [^1,^2]) and 
a random space direction, obtained by random triangulations of unit time strips [t, t + 1] 
by arbitrary but finite numbers of triangles with one edge along the time line t (resp. t + 1) 
and the opposite vertex on the time line t + 1 (resp. t). All other edges are then glued to 
their neighbors so as to form a triangulation. The boundary may be taken free, periodic 
or staircase- like on the left [H]. A typical such Lorentzian triangulation B in 1+lD reads 
as follows: 




> 

space 

These triangulations are best described in the dual picture by considering triangles as 
vertical half-edges and pairs of triangles that share a time-like (horizontal) edge as vertical 
edges between two consecutive time-slices. We may now concentrate on the transition 
between two consecutive time-slices which typically reads as follows: 



(2.1) 
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with say i half-edges on the bottom and j on the top (here for instance we have i = 9 
and j = 10). Denoting by \i) and i,j G Z the bottom and top states bases, we may 
describe the generation of a triangulation by the iterated action of a transfer operator 7 
with matrix elements Tjj- = (*^"'). Note that the corresponding matrix T = (Tij)ij^z+ is 
infinite. We shall deal with such matrices in the following. A useful way of thinking of 
them is via the double generating function: 



1 — z — w 

i,j>0 

2.2. Integrability. To make the model more realistic, we may include both area and 
curvature-dependent terms, by introducing Botlzmann weights w{Q) equal to the product 
of local weights of the form g per triangle (area term) and a per pair of consecutive triangles 
in a time-slice pointing in the same direction (both up or both down). The rules in the 
dual picture are as follows: 



g g a a 

For instance, in the example (12. ip above with i = 9 and j = 10, the product of local weights 
is g^^a^. It is easy to see that these weights correspond to a new transfer operator 7{g, a) 
with matrix elements 

Min(Mj) 

' i \ ( J \ -2k 



TM..,^ia,Y^' E (;)(y 



A:=0 

Equivalently, the double generating function reads: 



(2.2) fT(g,a){z,w) = T{g,a)ijz'w^ 



1 — ga{z + w) — g'^{l — o?)zw 

This model turns out to provide one of the simplest examples of quantum integrable 
system, with an infinite family of commuting transfer matrices. Indeed, we have: 

Theorem 2.1. [9] The transfer matrices T{g,a) and T{g',a') commute if and only if the 
parameters {g,a,g',a') are such that (f{g,a) = (f{g',a') where: 

l-g^{l-a') 
n9,a) = - 



ag 

This is easily proved by using the generating functions, and noting that for two infinite 
matrices A, B we have formally fAB = fA'*^fB, where * stands for the convolution product, 
namely 

(/a*/b)(^,w;) = jf ^^/a (^2;, fB{t,w) 
where the contour integral picks up the constant term in t. 
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This was extensively used[9j to diagonalize 7{g, a) and to compute correlation functions 
of boundaries in random Lorentzian triangulations. Our purpose here was simply to display 
the very simple form of the generating function fl2.2p . which will reappear later in these 
notes. 

2.3. Trees. For suitable choices of boundary conditions, the dual random Lorentzian tri- 
angulations introduced above may be viewed as random plane trees. This is easily realized 
by gluing all the bottom vertices of parallel vertical edges whose both top and bottom 
halves contribute to the curvature term (no interlacing with the neighboring time slices). 
A typical such example reads: 




Note that the tree is naturally rooted at its bottom vertex. 

To summarize, we have unearthed some integrable structure attached naturally to plane 
trees, one of the most fundamental objects of combinatorics. Note that in tree language 
the weights are respectively per edge, and a per pair of consecutive descendent edges 
and per pair of consecutive leaves at each vertex (from left to right). 

3. Planar maps and geodesics 

3.1. Two-dimensional quantum gravity. Two-dimensional quantum gravity (2DQG) 
is a theory describing the interactions of matter with the underlying space-time, in which 
both are quantized. Discrete models of 2DQG involve statistical matter models defined on 
statistical ensembles of discrete random surfaces, in the form of random tessellations. The 
Einstein action for these 2D random surfaces involves their two invariants, the area and the 
genus. The complete model will therefore combine the Boltzmann weights of the statistical 
matter model and of the underlying surface, in the form of parameters coupled to its area 
and its genus. Matrix integrals have proven to be extremely powerful tools for generating 
such discrete random surfaces with matter, by noting that the expansion for large matrix 
size matches the genus expansion. In a parallel way, the field theoretical descriptions of 
the (critical) continuum limit of 2DQG have blossomed into a more complete picture with 
identification of relevant operators and computation of their correlation functions. This was 
finally completed by an understanding in terms of the intersection theory of the moduli 
space of curves with punctures and fixed genus. Remarkably, in all these approaches a 
common integrable structure is always present. It takes the form of commuting flows in 
parameter space. 
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However, a number of issues were left unadressed by the matrix/field theoretical ap- 
proaches. What about the intrinsic geometry of the random surfaces? Correlators must be 
integrated w.r.t. the position of their insertions, leaving us only with topological invariants 
of the surfaces. But how to keep track say of the geodesic distances between two insertion 
points, while at the same time summing over all surface fluctuations? 

3.2. Maps and trees. Answers to these questions came from a better combinatorial un- 
derstanding of the structure of the (planar) tessellations involved in the discrete models. 
And, surprisingly, yet another form of integrability appeared. Following pioneering work 
of Schaeffer it was observed that all models of discrete 2DQG with a matrix model 
solution (at least in genus 0) could be expressed as statistical models of (decorated) trees, 
and moreover, the decorations allowed to keep track of geodesic distances between some 
faces of the tessellations. Marked planar tessellations of physicists are known as rooted 
planar maps in mathematics. They correspond to connected graphs (with vertices, edges, 
faces) embedded into the Riemann sphere. Such maps are usually represented on a plane 
with a distinguished face "at infinity", and a marked edge adjacent to that face. The degree 
of a vertex is the number of distinct half-edges adjacent to it, the degree of a face is the 
number of edges forming its boundary. 

Let us concentrate on the following example of tetravalent (degree 4) planar maps with 2 
univalent (degree 1) vertices, one of which is singled out as the root. The Schaeffer bijection 
associates to each of these a unique rooted tetravalent (with inner vertices of degree 4) tree 
called blossom-tree, with two types of leaves (black and white), and such that there is 
exactly one black leaf attached to each inner vertex. 




This is obtained by the following cutting algorithm: travel clockwise along the bordering 
edges of the face at infinity, starting from the root. For each traversed edge, cut it if and 
only if after the cut, the new graph remains connected, and replace the two newly formed 
half-edges by a black and a white leaf respectively in clockwise order. Once the loop is 
traveled, this has created a larger face at infinity. Repeat the procedure until the graph 
has only one face left: it is the desired blossom-tree, which we reroot at the other univalent 
vertex, while the original root is transformed into a white leaf. 
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Why is it a bijection? Remarkably, the information of black and white leaves is sufficient 
to close up the blossom-tree into a unique map. Clockwise around the tree, match each 
black leaf to its immediate follower if it is white, and repeat until all but one leaves are 
matched. The loser of the musical chairs is the root of the map. 



3.3. Exact enumeration and integrability. The bijection above allows to keep track 
of the geodesic distance between the faces containing the two univalent vertices, defined as 
the minimal number of edges crossed in paths between the two faces (this distance is 2 in 
the example). Let us introduce a weighted counting of the maps, by including a weight g 
per vertex. 

Defining Rn{g) to be the generating function for maps with geodesic distance < n between 
the two univalent vertices, we have the following relation: 

(3.1) Rn{9) = 1 + gRn{9) {Rn+M + K(^?) + Rn-M) 

easily derived by inspecting the environment of the vertex attached to the root of the tree 
when it exists. 

This recursion relation must be supplemented with boundary conditions. First, when 
n = 0, the relation (13. ip makes sense only if the term R^i{g) is omitted, so we set R-i{g) = 
0. Moreover, if we take n — >■ oo, we simply relax the distance condition, and the function 
R{g) = limn^^Rn{g) is the generating function for maps with two univalent vertices. It 
satisfies the limiting equation R = 1 + 3gR^, and we easily get 

(3.2) R(,) ^ l^^l^ 

6g 

as the unique solution with a formal power series expansion of the form R{g) = 1 + 0{g). 

The function Rn{g) is the unique solution to (13. ip for n > such that (i) R-i{g) = 
and (n) limn^oo Rn{g) = Rid) of (13. 2p . A first remark is in order: the equation (13. ip . 
viewed as governing the evolution of the quantity Rn{g) in the discrete time variable n, 
is a classical discrete integrable system. By this we mean that it has a discrete integral of 
motion, expressed as follows. The function (j){x, y) defined by 

(3.3) (t){x,y) = xy{l-g{x + y))-x-y 

is such that for any solution Sn of the recursion relation (13.11) . the quantity 0(S'„,, Sn+i) is 
independent of n. In other words, the quantity 0(5'„, Sn+i) is conserved modulo (13. ip . This 
is easily shown by factoring </)(S'„, Sn+i) — 4>{Sn-i, Sn)- 
By the conservation of (f){Rn{g) , Rn+i{g)) , we find that 

(j){Rn{g),Rn+i{g))= lim (f){RU9),Rm+i{9)) = <I>{R{9),R{9)) 

which gives us an explicit relation between Rn{g) and Rn+i{g)- It turns out that we can 
solve explicitly for Rn{g)'- 
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Theorem 3.1. [8J The generating function Rn{g) for rooted tetravalent planar maps with 
two univalent vertices at geodesic distance at most n from each other reads: 

o _ p( A^-^i9r^')i^-x{gr') 

where x{g) is the unique solution of the equation: 

1 1 

x-\ h4 



X gR{.gf 

with a power series expansion of the form x{g) = g + 0{g^). 

The form of the solution in Theorem 13.11 is that of a discrete sohton with tau-function 
Tn = I — x{gY. Imposing more general boundary conditions on the equation (13.11) leads 
to elliptic solutions of the same flavor. The solution above and its generalizations to many 
classes of planar maps have allowed for a better understanding of the critical behavior 
of surfaces and their intrinsic geometry. Recent developments include planar three-point 
correlations, as well as higher genus results. 

To summarize, we have seen yet another integrable structure emerge in relation to (dec- 
orated) trees. This is of a completely different nature from the one discussed in Section 
[21 where a quantum integrable structure was attached to rooted planar trees. Here we 
have a discrete classical integrable system, with soliton-like solutions. Such structures will 
reappear in the following, in relation to (fixed) lattice statistical models. 

4. Alternating sign matrices 

4.1. Lambda-determinant and Alternating Sign Matrices. The definition of the so- 
called Lambda-determinant of Robbins and Rumsey [32] is based on the famous Dodgson 
condensation algorithm [17] for computing determinants, itself based on the Desnanot- 
Jacobi equation, a particular Pliicker relation, relating minors of any square k + 1 ^ k + 1 
matrix M: 

(4.1) |M| X \Ml'',ll\ = \M',+l\ X \Ml\ - |M^l| x \Ml^,\ 

where \M^^f^'"'fj'\ stands for the determinant of the matrix obtained from M by deleting 
rows ii, ...,ir and columns ji, The relation (14.11) may be used as a recursion relation 

on the size of the matrix, allowing for efficiently compute its determinant. 

More formally, we may recast the algorithm using the so-called A^o T-system (also known 
as discrete Hirota) relation: 

(4.2) Tij^k+iTij^k-i = Tij^i^kTij-i± — Tj-|_ij^fcTj_ij^fc 

for any i,j, k E Z with fixed parity of i + j + k. Now let A = (ajj)ijg{i^2,...,n} be a fixed 
n X n matrix. Together with the initial data: 

Te,m,o = 1 {£, m E Tj] i + m = n mod 2) 

(4.3) Tij^i = Q j-i+n+i i+j+n+i {i,jEZ]i+j = n + lm.od2]\i\ + \j\<n — l), 



2 ' 2 
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the solution of the T-system (14.21) satisfies: 

(4.4) To,o,„ = det{A) 

Given a fixed formal parameter A, the Lambda-determinant of the matrix A, denoted by 
\A\x is simply defined as the solution To^o,n = |^|a of the deformed T-system 

(4.5) Tij^k+iTij^k-i = Tij^i^kTij-i^k + A Tj+i j'^fcTj_i j^fc 

subject to the initial condition (14. 3p . 

The discovery of Robbins and Rumsey is that the Lambda-determinant is a homogeneous 
Laurent polynomial of the matrix entries of degree n, and that moreover the monomials in 
the expression are coded hj n x n matrices B with entries bij G {0, 1, —1}, characterized 
by the fact that their row and column sums are 1 and that the partial row and column 
sums are non-negative, namely 

k k 

>0 Y.^J,^>0 ik = l,2,...,n-l;j = l,2,...,n) 

i=l 1=1 
n n 

5^6,, = 1 J2^j, = l (j = l,2,...,n) 

i=l i=l 

Such matrices B are called alternating sign matrices (ASM). These include the permutation 
matrices (the ASMs with no —1 entry). Here are the 7 ASMs of size 3: 

^0 1 
1 -1 
.0 1 

There is an explicit formula for the Lambda-determinant [35] : 









(4.6) \AU= Yl A^'^^(^)-^(^)(l + A)^(^) n«'T' 



'■J 

nxnASMB i,j 



where Inv(i?) and N{B) denote respectively the inversion number and the number of entries 
— 1 in S, with 



Inv(5) = Y Kibj,k 



1 < i < J < n 
l<fc<^<ri 



N{B) = -{-n+ J2 

l<i,j<n 



Note that for A = —1, only the ASMs with N{B) = contribute, i.e. the permutation 
matrices, for which Inv(i?) coincides with the usual inversion number of the corresponding 
permutation, and therefore (14. 6 p reduces to the usual formula for the determinant. 

Robbins and Rumsey formulated the famous ASM conjecture, that the total number of 
nxn ASMs is given by A„ = JJpo fw ■ 
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Figure 1. The combinatorial family of ASMs. From left to right: ASM, 6V- 
DWBC and FPL , all in bijection; dense loop gas: its groundstate/limiting proba- 
bility vector satisfies the qKZ equation, the components measure FPL correlations 
(RS conjecture), their sum matches the 6V-DWBC partition function with in- 
homogeneous spectral parameters Zi and = 1; DPP: their refined evaluation 
matches that of ASMs (MRR conjecture); TSSCPP: their refined enumeration 
matches a sum rule for qKZ solutions at generic q and = 1; Variety = 0: its 
degree/multidegree matches solutions of qKZ for q = 1. 



4.2. Integrabilities: six vertex, loop gas and more. ASMs relate to two different 
kinds of integrable systems. 

On one hand, the T-system relation (14. 2 p and its deformations (14. 5 p are known to be dis- 
crete classical integrable systems [27], when interpreted as a 2+lD evolution in the discrete 
time variable k. It is also part of a combinatorial structure called a cluster algebra [20], 
which we will discuss in the next section. As such, it has a guaranteed Laurent property, 
namely that its solutions are always Laurent polynomials of the initial data, which explains 
the miracle witnessed by Robbins and Rumsey. 

On the other hand ASMs are in bijection with configurations of a famous integrable 2D 
lattice model, the Six Vertex (6V) model with particular, so-called Domain Wall Boundary 
Conditions (DWBC) |25l|23]. This latter connection allowed Kuperberg to derive an elegant 
proof of the ASM conjecture [2S], shortly after Zeilberger's remarkable direct proof using 
generating functions [39j. 

The connection with the 6V model is only the tip of an iceberg, which has kept revealing 
more structure over the years (see the sketch in FigH]). The ASMs of size n are indeed also 
in same number as: 

(i) the Totally Symmetric Self- Complementary Plane Partitions (TSSCPP), namely the 
rhombus tilings of a regular hexagon of size 2n drawn on the triangular lattice, with all 
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possible symmetries of the hexagon, which amount to the tihngs of a fundamental domain 
occupying ^ of the hexagon (see FigH]). 

(ii) the Descending Plane Partitions (DPP) of order n, which we will define below. 

There is however no known natural bijection between these objects, and it remains a 
challenge to find one. 

The picture has grown more complex with the remark of Batchelor, DeGier and Nienhuis 
that the groundstate of the so-called 0{n = 1) dense loop model on a cylinder of perimeter 
2n has positive rational entries summing up to 1 with common denominator An. Alter- 
natively, each component is interpreted as the limiting probability of loop configurations 
realizing a fixed pattern of connection of the 2n string ends on the boundary of the cylinder. 
Shortly after, the numerator of each component of this vector has been identified by Razu- 
mov and Stroganov (RS) as some particular correlation of the 6V model (and equivalently 
of ASMs), once reformulated as a model of fully-packed loops (FPL) on a n x n square 
lattice grid [31]. This conjecture was proved recently by Cantini and Sportiello [7j in a 
purely combinatorial way. 

In the next subsections, we will address the connection between refined enumeration 
of ASMs and DPPs (the Mills, Robbins and Rumsey (MRR) conjecture [33]). Before 
we go into this, let us mention that more connections have been found between the var- 
ious combinatorial objects. The computation of the groundstate of the dense loop gas 
with multi-parameter inhomogeneous Boltzmann weights was performed [H] by use of the 
so-called quantum Knizhnik-Zamolodchikov equation, whose suitably normalized solution 
yields the polynomial components of the groundstate vector, depending on 2n complex 
spectral parameters Zi, i = l,2,...,2n and one quantum parameter q parameterizing the 
weight per loop —{q + q~^)- it was shown [H] that the sum of these components equates 
the partition function of the 6V model with DWBC when g is a non-trivial cubic root of 
unity, expressed as a determinant by Izerginp3]. The ASM number is recovered in the 
limit when Zi ^ 1 for all i. On the other hand, when all Zi = 1 but q arbitrary, it was 
shown that the sum of the components of the groundstate equates the partition function 
for TSSCPP on the fundamental domain above, with a particular weight {q + q~^) for one 
type of tile [TB] . The solutions of the qKZ equation therefore form a bridge between ASM 
and TSSCPP, in that two different specializations relate to either objects. 

Finally, last but not least, yet another combinatorial problem turns out to be related 
to all of the above. It was shown in [16] that the degree of the variety of strictly upper 
triangular 2n x 2n complex matrices of vanishing square is equal to this weighted TSSCPP 
counting function, for q = 1. Moreover, the variety actually decomposes into irreducible 
components whose degrees are the entries of the vector, solution to the qKZ equation in the 
correlated limit where g — )■ 1 and Zi ^ 1. More generally this extends to the multi-degree 
of this variety, by considering its torus-equivariant cohomology, and by taking a suitable 
limit of the qKZ solution [TBI ES] • 

4.3. ASMs from 6V. The 6V-DWBC model on a square n x n grid has configurations 
obtained by choosing an orientation for each edge of the square lattice, in such a way 
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that the flow is conserved at each vertex, namely exactly two edges are incoming and two 
outgoing at each vertex. Moreover the Domain Wall Boundary Conditions impose that 
horizontal external edges to the grid point into the grid, while vertical ones point out of 
the grid. This gives rise to the following 6 possible local configurations: 



6V 








h + 


6V 
DWBC 


ASM ( 


[) 


1 -1 






where we also have indicated the corresponding ASM entries. One usually attaches Boltz- 
mann weights to each of the above configurations as follows: each row (resp. column) of 
the grid carries a complex number Zi, i = 1,2, ...n (resp. Wi, i = 1,2, ...n ) called spectral 
parameter. Moreover the weights depend on a "quantum" parameter g G C*. We have the 
following parametrization of the weights: 



a{z,w) = qz — q w b{z,w) = q z — qw c{z,w) = {q —q )\fzw 

where a(^z,vS) is the weight for a vertex of type a\ or at the intersection of a line with 
parameter z and column with parameter w, etc. With this parameterization, the model has 
an infinite family of commuting row-to-row transfer matrices, and can be exactly solved by 
Bethe Ansatz techniques. Using recursion relations of Korepin[25j, lzergin[23j obtained a 
compact determinantal formula for the partition function of this 6V-DWBC model, defined 
as the sum over edge configurations of the product of local vertex weights, divided by the 
normalization factor nr=i Wj) (to make the answer polynomial in the z's and w's). It 
reads: 

An), r . r .^ Ui,, (^{Zi, Wj)b{Zi, Wj) 



(4.7) {--.},{.,}) - -%(,)^(,) - .det^. 

where A{z) = Y[i<i<j<ni^i ~ stands for the Vandermonde determinant of the z^s. 

In the above bijection between 6V configurations and to ASMs, it is easy to track both 
quantities N[B) and lnv{B) in terms of 6V weights. We find that 



N - n N 
N{B) = N,, = Inv(5) - N{B) = N^, = N,, = ^ 



where Na-, N^-, N^^ stand for the total numbers vertex configurations of each type. The 
determinant result above can therefore be used to compute the refined partition func- 
tion for ASMs, which counts ASMs with a weight x/y per entry —1 and a weight y for 
each inversion (with these weights, the 7 ASMs of size 3 above have respective weights: 

Setting 

/C\ 2 /ON 2 

(4.8) X -- 



bJ " \b 
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we have: 

(4.9) 41,(x, y)= Y: = b, c) 

ASM B 

where Zgy (a, b, c) refers to the homogeneous hmit of the partition function (14. 7p of the 6 V 
model in which all a{zi, Wj) are equal to a, etc. (This is obtained by letting all Zi ^ z and 
all Wi — )■ w, with a = a{z,w), b = b{z,w) and c = c{z,w).) This and more refinements 
were worked out in [3]. We have the following remarkable result: 

Theorem 4.1. |3] The partition function for refined ASM reads: 

(4.10) 4"J^(x,i/)= det (il-u)I + uG) 

0<i,j<n— 1 

where v is any solution to the equation 

(4.11) xiy{l - u) = u + y{l - u) 

and the n x n determinant is the principal minor for the n first rows and columns of the 
infinite matrix Masm = (1 ~ '^)I + '^G whose entries are generated by 

1 — V V 
1 — zw 1 — zx — w — [y — xjzw 

We notice that the second term in the generating function is nothing but that of the 
transfer matrix for 1+lD Lorentzian gravity, up to some gauge transformation z — )■ zj^fx 
and w — )■ w^fx^ and upon the identification x = a^g"^ and y = g^. 

4.4. DPP from lattice paths. Descending plane partitions are arrays of positive integers 
of the form: 

0-1,1 0'1,2 Ctl,3 ■■■ ■■■ 0^1,Ati-2 

^2,2 ^2,3 ■ ■ ■ ■ ■ ■ '^2,fj.2 

(^3,3 ■ ■ ■ "^S./is 



(^r,r ' ' ' Or,^r 

such that the sequece /ij is strictly decreasing /ij+i < /ij, and that for Aj = /Uj — ? + 1, 
Ao = oo: 

^ (^i,j+l Oi,j > Oi+l,j -^i ^ Oi,i — ^i-1 

for all i,j. 

The integers j are called parts. A DPP is said to be of order n if aj j < n for all i,j. A 
part j is said to be special if a^j- < j — i. By convention, the empty partition is a DPP. 
In the following we will compute the partition function Z^j^pp{x,y) which is the sum over 
all DPP of order n with a weight x per special part and y per remaining part. Here are 
the 7 DPP of order 3: 

2 3 31 32 33 ^ ^ 
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1 2 3 4 5 6 



Figure 2. Non-intersecting lattice path configuration for a sample DPP of order 
n > 8. We have indicated the domains in which horizontal steps correspond to 
special/non-special parts. 



They have respective weights: l,y,y, xy, y'^, y"^, y^, as the only special part is the entry 1 in 
the fourth DPP. 

The DPPs are in bijection with configurations of non-intersecting lattice paths illustrated 
in Figj2] and defined as follows. We consider the non-negative quadrant Z_|_ x Z_|_ of the 
integer plane. Paths may take two types of step: vertical up or horizontal left. They start 
along the x axis at positions of the form (sj, 0) (recorded from right to left) and end along 
the y axis at positions (0, Si + 2) (recorded from top to bottom) for i = 1, 2, r. We add 
a final horizontal left step at the end of each path. Reading paths from left to right and 
top to bottom, we record the vertical positions y = aij of the j-th horizontal step from 
the left taken on the i-th path from top (steps with y = are not recorded). These form a 
DPP with r rows, of order any n > Si + 2. Conversely to each DPP with r rows we may 
associate such a path configuration. Note that the starting points are such that Si = Aj — 1, 
where Xi = fii — i + 1 the total number of parts in the row i. We note that the special parts 
correspond to horizontal steps taken in the strict upper octant y > x + 1 of the plane, and 
the remaining parts correspond to the horizontal steps in the domain I < y < x + 1, while 
horizontal steps along the x axis do not count. 

The computation of Z^^pp{x^y) uses the Lindstrom-Gessel-Viennot [321 [22] determinant 
formula expressing the partition function for non-intersecting lattice paths Z as a determi- 
nant det(Zjj) where Zi^j is the partition function for a single path starting from the i-th 
point and ending at the j-th. This leads to the following: 

Theorem 4.2. [3j The partition function Z^^pp{x, y) for DPP of order n with weight x per 
special part and y per other part reads: 

(4.13) Z^j^lp{x,y)= det (I + H) 
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where the determinant is the principal minor of size n x n of the infinite matrix M^pp = 
I + H with generating function: 

(4.14) /aw(^, w) = + ^ — 

1 — zw 1 — zl — xz — w — [y — xjzw 

This bears a striking resemblance to the formula of Theorem 14.11 and also involves the 
generating function for the transfer matrix of 1+lD Lorentzian gravity. 

4.5. Proof of the ASM-DPP conjecture. The complete MRR conjecture is now proved[3]. 
The proof uses the same ingredients as exposed here, and includes extra refinements. For 
pedagogical reasons, we retain here only the special/non-special part dependence on the 
DPP side, and the number of — 1/inversion number dependence on the ASM side. Let us 
show that: 

(4-15) Z^ASMi.^^ y) = Zolpix, y) 

Both expressions are determinants of the principal minor of size n of some infinite matrix, 
in other words, these are the determinants of a finite truncation to the n first rows and 
columns of infinite matrices. There is a very simple relation (independent of n) between 
the generating functions of the two infinite matrices Masm and M^pp, namely: 

(1 - Z)(l - (1 - v)w)fMaPp{z,w) = (1 - - w)fMASM{z,w) 

clS cl direct consequence of (14. lip . 

Moreover, introducing the strictly lower triangular infinite "shift" matrix S with entries 
Si^j = Sij+i for i,j > 0, we see that for any infinite matrix A and any complex parameters 
a,/3: (1 - az)fA{z,w) = f{i-aS)A{z,w), and similarly (1 - (3w)fA{z,w) = fA{i-f}s'){z,w). 
Moreover, as I — aS is lower uni-triangular and I — f3S^ is upper uni-triangular, any finite 
truncation to the n first rows and columns of (I — aS)A and A{I — fiS^) has the same de- 
terminant as the corresponding truncation of A alone. This allows to identify the principal 
minors of M^pp and Masm and the desired identity follows. 

To conclude, our generating function method has allowed for a simple comparison and 
eventual proof of the identity of the partition functions for ASM and DPP. In passing, 
it shows some mysterious connection to the quantum integrable transfer matrix of 1+lD 
Lorentzian gravity. It seems that the quantum integrability of the 6V model is unrelated to 
this one, and this is yet another puzzling fact about ASMs as a crossroad between various 
integrable models. 

5. Discrete Integrable systems and cluster algebras 

5.1. T-system and initial data. Let us concentrate on the T-system realtion (14. 5 p with 
for simplicity A = 1. We may use this relation to describe the non-linear discrete time 
evolution (/c G Z) of a 2D space-dependent quantity (for i,j G Z), Tj j fc, with: 

l-^-J-J '^i,j,k+l = Tf: '■ 

i,j,k—l 
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Figure 3. The "flat" stepped surface carries all the initial data at times k = 
0, 1 for the T-system. We have represented a local application of the T-system 
evolution: it evolves the stepped surface by "adding" an elementary octahedron 
to it. 

Initial conditions determining Tij^k entirely consist in assigning values of Tij^k along any 
fixed "stepped surface" k = kij)}ij^z with kij G Z such that \kij+i — kij\ = — 
kij\ = 1 for all i,j G Z. The simplest such surface, denoted by kg, has kij = i + j mod 2, 
and consists of all the points at times and 1 (cf. FigJS]). This is the surface used in the 
definition of the Lambda-determinant. The evolution (15.11) may be used to trade a stepped 
surface k with a surface k' identical to it except for the point {i,j,kij) — )■ {i,j,ki j = 
kij + 2). This is however possible only if the four neighboring points share the same value 
of A; = ki_i j = fcj+ij- = kij+i = kij_i = kij + 1. We see that this amounts to completing 
an octahedron and trading its point with the lowest value of k for the one with the largest 
(see Fig|2])- The application of the T-system is sometimes called "octahedron move". 

This will be later interpreted as a mutation at point in the corresponding cluster 
algebra. 

5.2. Cluster Algebra. Here we give a simplified definition of cluster algebra (namely 
skew-symmetric or geometric type, with no coefficient), and we refer to Fomin and Zelevin- 
sky's original paper^20| for more general considerations. 

Let T be an infinite tree with fixed vertex degree n, and edges labeled 1,2, ...,n around 
each vertex. To each vertex f of T we attach two data: (i) a vector of formal variables with 
n components = (xi(f), ...,x„(f)) (ii) an integer skew-symmetric matrix = (6jj(f)), 
coding a quiver (graph with oriented edges) with n vertices labeled 1, 2, n, and with 
no oriented loops of length < 2. 

The data at adjacent vertices v,v' of 7 connected via an edge labeled k are related via 
a mutation /x^. The mutation acts both on quivers and data vectors. 
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The application of fik on transforms it as follows: (i) it reflects all edges of adjacent 
to the vertex k (ii) for any path of length 2 of the form z — )■ A; — )■ j in Q^,, it adds an oriented 
edge i — )■ j to Qp^^{v) unless some edges j ^ i already exist, in which case one of them is 
canceled. As an example, we represent here the mutations fii and then /i2 on a quiver Q: 

G) K4) G> (4) (i> ® 





^2 




The application of fi^ on leaves Xj 
(5.2) Xkiv)xk{v') = 



XAV 



k, where 



n 



i— >fc in ( 



invariant, except for i 

k—>-j in Qv 

In other words, the quiver at v encodes the transformation of the data vector from 
f — f ' any adjacent vertex on T. 

The cluster algebra is the algebra with a unit generated by all the cluster variables at 
all the clusters. It is entirely determined by the "initial data" at any fixed vertex of the 
tree 7. The data vectors are called clusters, their entries cluster variables. The matrix B 
is called exchange matrix. The pair (x, B) is called a seed. The integer n is called the rank 
of the cluster algebra, and can be infinite. 

5.3. Properties, applications, conjectures. The major built-in property of cluster alge- 
bras is the Laurent phenomenon: the cluster variables at any vertex of T may be expressed 
as Laurent polynomials of the cluster variables at any other vertex of the tree. This is not 
obvious, as the mutations involve many divisions. The fact that polynomials are always 
divisible by the new denominators is a non-trivial outcome of the definitions. 

Two important results have been obtained regarding the classification of cluster alge- 
bras. First the cluster algebras with finitely many distinct cluster variables (finite type) 
are classified [21]. In our case, they correspond to an initial quiver Q that is the Dynkin 
diagram of a simply-laced Lie algebra of A, D, E type, with a choice of orientation of the 
edges. Next, the cluster algebras with finitely many distinct quivers (mutation-finite) are 
also classified [19]. They are essentially made of cluster algebras attached to (decorated) 
triangulations of Teichmiiller space, plus a few exceptions. 

Applications of cluster algebra are very diverse both in mathematics and physics: canon- 
ical and dual canonical bases of quantum groups, total positivity, preprojective algebra 
representations, quiver representations, triangulated 2-Calabi Yau categories, discrete in- 
tegrable systems and somos-like sequences, dimer models, Donaldson-Thomas invariants of 
topological string theory, wall-crossing in super- Yang-Mills theories, quantum dilogarithm, 
Teichmiiller space geometry, etc. 

A mysterious property was conjectured by Fomin and Zelevinsky[20], and was only proved 
in particular cases so far: the Laurent phenomenon of cluster algebra is positive, namely 
the Laurent polynomials involved have only non-negative integer coefficients. This calls for 
some combinatorial interpretation, still missing to this day. 
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From now on, we concentrate on the connection to discrete integrable systems, and more 
precisely to the T-system. 

5.4. T-system as a sub-cluster algebra. We have the following: 

Theorem 5.1. [lOj The T-system relation is a mutation in an infinite rank cluster algebra. 

The initial seed is constructed as follows. The cluster is made of the initial data along 
the stepped surface ko discussed above, namely with x = {xij)ij^z, and Xij = Tjj_j+jmod2- 
The quiver may be represented as an infinite square lattice with vertices {i,j) G and 
oriented edges (i, j) — )■ {i,j ± 1) if i + j = mod 2, and {i,j) {i ± otherwise. 

Starting from the initial data stepped surface kg, we may apply a mutation /ij whenever 
the octahedron move is permitted at (i, j, kij) (this limits the number of possible mutations, 
and therefore we do not explore the entire cluster algebra in this way). It is easy to check 
that in all cases the quiver has 2 incoming and two outgoing edges at each mutable vertex, 
and the rules 05.21) are equivalent to the T-system relation (15.11) . (This can be inferred from 
the example of mutation given above.). 

5.5. Solution and Laurent positivity. Let us briefly describe the solution of the T- 
system in the case of the initial data surface kg. The surface can be decomposed into 
half-octahedra that can be completed upon mutation. We color in white and grey the 
triangular faces of each half-octahedron as follows, and attach to the "rhombi" formed by 
pairs of triangles sharing a horizontal edge (belonging to the [j, k) plane) the following 
2x2 matrices D, U according to whether the grey triangle is on the bottom or the top of 
the rhombus: 



V 




u 



where a = Tij_i^k, b = Tij^k-i, c = Ti^j+i,k u = Ti^ij^k and v = Tj+ij-fc and where the 
subscript i indicates that we should think of these 2x2 matrices as embedded into an 
infinite identity matrix with rows and columns indexed by Z, except at rows and columns 
i,i + 1 where we insert the 2x2 matrix. 

The octahedron move corresponds to the matrix identity 




= Ui{a, X, v)Di{u, X, c) 



u 
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which is satisfied ifi^ 

bx = ac + uv 

By this construction, we represent the T-system mutation by a reordering relation of D 
and U matrices. 

The solution Tij^k for A; > and i+j + k = Oof the T-system with inital data along ko 
can be expressed explicitly in terms of the assigned values as the determinant of a matrix 
Mj J built out of products of D's and t/'s. The product is best expressed pictorially: 




= A^/^A-l^/i-lA+l^/^+l...A^/^ 



where the picture is a view from behind of the stepped surface ko ( from a point k') 
with k' « 0), the central vertex is (i, j,i+j mod 2) and the size of the tilted square is k — 1, 
so that all the boundary vertices sit at time k = 1. Moreover, the product is read from 
left to right, and a D,U is to the left of a D,U in the product iff it is on the picture, with 
the obvious fact that A, Ui commute with Dj, Uj as soon as |^ — j | > 1. Note that Mij^k 
is the embedding of a 2A; — 2 x 2A; — 2 matrix into an infinite identity matrix. Truncating 
{Mij^k)a,b to the square window a,b E [i — k + 2,i + k — 1] we may take its determinant, 
still denoted det(Mjj jt). 
We have the final: 

Theorem 5.2. [13j The solution Tij^^ of the T-system is expressed in terms of the initial 
data along the stepped surface ko as: 

k-l 

TiJ,k = det(Mij-fc) Y^'^i'-k+l+aJ-a^Ti-k+l+aj+a,! 
a=l 

This provides an explicitly positive Laurent polynomial of the initial data, due to the 
structure of the D, U matrices. The proof is by induction under octahedron move, making 
extensive use of the relation 05.41] . This was extended to arbitrary stepped surface initial 
data[13j as well as other boundary conditions such as wall i =const. or j =const. along 
which Tj J- is fixed to be 1. In all cases positivity follows from explicit expressions involving 
the D, U matrices. 

6. Conclusion 

In these notes we have explored various combinatorial problems, with one or more links 
to integrable systems. These links seem to always point to the existence of some kind of 
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tree or path model underlying the original combinatorial problem. Note that a path model 
is typically a discrete version of a quantum path integral. 

In the last section, the solution of the T-system was expressed in terms of U matri- 
ces. These can be interpreted as elementary pieces allowing to construct networks, namely 
directed graphs with edge weights. More precisely, we attach to Di, Ui the following ele- 
mentary pieces of graph, connecting entry and exit vertices labeled + l via edges oriented 
from left to right, and whose weights match the corresponding matrix entries: 



^ it i+1 



Di{u,a,b) 




Ui{u,a,b) 




The above products Mjj ^ of D,U matrices correspond to the following networks: 

i+3 » >> ■ .> 




with a total of 2k — 2 rows. By the Lindstrom-Gessel-Viennot theorem, the factor det{Mij^k) 
in Theorem 15. 2 1 is the partition function for k non- intersecting paths starting at the k south- 
west entry points of the network, and ending at the k south-east exit points of the network. 
This gives yet another path model for the solution of the T-system. The positivity of the 
path weights is responsible for the positive Laurent phenomenon. 

Path models are also very important for non-commutative generalizations. Indeed, a 
path is a fundamentally non-commuting object, as its steps are taken in a particular order. 
Attaching non-commuting weights to paths yields natural non-commutative expressions. 
Cluster algebras have a "quantum" version, in which the cluster variables have commu- 
tation relations encoded by the quiver [5]. In rank 2, a fully non-commutative version 
of cluster algebra was introduced via path model solutions pTl [T2] . Finally, a higher rank 
non-commutative discrete Hirota equation was derived [12] from non-commutative path and 
continued fraction expressions involving quasi-determinants of non-commuting variables. 
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